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with  a  global  model  which  Is  based  on  the  AFGL  baseline  model.  A  set  of  eleven  96  hour 
forecast  were  run.  In  each  case  two  experiments  were  conducted,  one  in  which  the  sea 
surface  temperatures  were  fixed  at  their  climatological  monthly  mean  values  and  one  in 
which  the  temperature  of  the  eastern  Mediterranean  was  Increased  by  up  to  2.5  degrees. 

In  the  control  runs,  it  was  found  that  the  rainfall  in  this  region  was  produced  primarily 
by  the  convective  scheme  of  themodel.  In  the  anomaly  runs,  the  rainfall  incresed  mainly 
over  the  area  of  the  imposed  anomaly  and  over  the  surronding  region.  The  increase  in 
ralnfalliwas  traced  to  the  combined  effects  of  additional  moisture  availability  due 
to  enhaced  evaporation  and  increasediconvective  activity  due  to  the  destabilizing  effect 
of  heating  from  below  by  the  enhanced  sensible  heat  flux. 

A  detailed  study  was  made  of  the  6  November  1986  case  in  which  we  ran  several  experiments 
with  different  anomaly  magnitudes  and  with  different  model  resolutions.  The  overall 
agreement  between  the  low  and  high  resolution  runs  la  very  encouraging.  There  are  of 
course  differences  in  the  details  due  to  the  ability  of  the  higher  resolution  model  to 
better  represent  smaller  scale  features  in  the  humidity  field. 
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1.  INTRODUCTION 


Air  sea  interaction  over  the  eastern  Mediterranean  plays 
an  important  role  in  shaping  the  weather  of  the  Middle  East. 
In  this  region,  the  year  can  be  divided  into  two  main 
seasons,  a  hot  dry  summer  from  April  to  September  and  a  cool 
wet  winter  from  October  to  March.  Precipitation  in  this  area 
is  usually  associated  with  migratory  cyclones  which  cross 
the  Mediterranean  and  bring  cold  air  masses  from  Europe 
(Tzvetkov  et  al.,  1985).  On  occasion,  a  secondary  low  will 
form  over  the  sea  in  the  vicinity  of  Cyprus.  These  so  called 
Cyprus  lows  will  often  enhace  the  precipitation  in  this 
region.  Ozsoy  (1981)  noted  that  cyclogenetic  activity  in 
this  region  is  also  associated  with  strong  interaction 
between  the  Subtropical  and  the  Polar  Front  Jet  streams. 

The  seasonal  rainfall  over  the  eastern  Mediterranean 
is  heaviest  along  the  northern  and  eastern  coasts  and 
generally  decreases  southward  and  eastward.  In  a  statistical 
study  of  the  region,  Tzvetkov  and  Assaf  (1982)  found  a 
remarkably  high  correlation  of  0.9  between  winter 
precipitation  in  Israel  and  what  they  termed  heat 
storage  extraction  from  the  eastern  Mediterranean.  While 
these  results  clearly  demonstrate  the  Importance  of  the 
sea  they  do  not  explain  the  physical  mechanism  responsible. 

In  this  report,  we  have  taken  the  first  step  in 
understanding  the  complex  physical  processes  that  control 
rainfall  In  this  region.  In  order  to  accomplish  this, 
we  have  used  a  modified  version  of  the  global  model 
developed  at  the  Air  Force  Geophysics  Laboratory  (AFGL).  In 
Section  2  we  describe  the  changes  that  were  implemented  in 
the  model  while  in  Section  3  we  present  results  from  a 
series  of  experiments  in  which  the  sea  surface 
temperatures  in  the  eastern  Mediterranean  were  varied  from 
their  climatological  mean  values.  Section  4  contains  our 
cone  1  us  1 ons . 


2.  MODIFICATIONS  TO  THE  AFGL  BASELINE  MODEL 


The  AFGL  baseline  model  has  already  been  described  in 
detail  by  Brenner  et  al.  (1982,  1964)  and  will  only 
be  briefly  summarized  here.  It  is  a  global  primitive 
equation  model  on  a  longitude,  latitude,  sigma  grid.  The 
horizontal  representation  is  handled  spectrally  while  the 
vertical  differencing  Is  based  on  the  quadratic 
conserving  scheme  described  by  Arakawa  and  Lamb  (1977). 
The  time  scheme  is  semi - imp  1 Icl t .  Subgrid  scale  physical 
parameter  1 zatlons  are  based  on  the  schemes  used  in  the 
previous  version  of  the  National  Meteorological  Center 
(NMC)  global  spectral  model  (Sela,  1980).  These  include: 
(1)  bulk  aerodynamic  formulae  for  surface  drag,  sensible 
heat  flux,  and  evaporation  (the  latter  two  over  the  oceans 
only))  (11)  Kuo  (1965)  type  moist  convection)  (Ill)  large 
scale  condensation;  (iv)  dry  convective  adjustment)  and  (v) 
horizontal  diffusion.  In  this  section  we  will  describe  the 
modifications  that  we  made  to  the  baseline  model  for  the 
purposes  of  our  study.  Our  efforts  focused  mainly  on 
improving  the  planetary  boundary  layer  scheme  of  the  model. 


2.1  Surface  fixed  fields 


The  surface  fixed  fields  in  a  numerical  model  describe 
the  physical  structure  of  the  surface  of  the  earth.  These 
form  an  integral  part  of  any  scheme  designed  to  simulate  the 
interaction  between  the  atmosphere  and  the  solid  earth.  For 
the  AFGL  baseline  model,  these  fields  include  topography, 
land  sea  mask,  sea  surface  temperature,  and  the 
neutral  drag  coefficient.  A  procedure  for  deriving 
consistent  surface  fixed  fields  for  this  model  has  been 
described  by  Brenner  (1987)  and  will  only  be  briefly 
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reviewed  here.  The  method  is  based  on  the  assumption  that  the 
basic  data  Is  available  on  a  2.5*  X  2.5*  latitude 
longitude  grid.  The  required  values  for  the  model  are 
derived  from  this  basic  data  set  by  Interpolating  onto 
the  model  grid  with  a  procedure  that  Is  similar  to  one 
pass  of  a  Cressman  (1959)  type  of  objective  analysis. 

The  AFGL  baseline  model,  as  many  other  models,  uses  a 
smoothed  topography  which  Is  generated  by  subjecting  the 
2.5*  X  2.5*  data  to  various  smoothing  operators  and  then 
Interpolating  onto  the  model  grid.  Recently  It  has  been 
recognized  that  using  a  rougher  topography  could  be 
beneficial  In  Improving  medium  range  forecasts.  The  attempt 
to  Incorporate  the  effects  of  the  subgrid  scale  variance  of 
the  topography  has  resulted  in  various  concepts  such  as 
envelope  orography  (Wallace  et  al.,  1983)  and  step 
mountains  (Meslnger,  1985).  However,  the  "best" 
representation  of  topography,  especially  In  spectral 
models,  has  not  yet  been  established.  For  our  study  we 
have  adopted  an  approach  which  Is  somewhere  between  the 
smoothed  topography  and  the  various  rough  topographies 
mentioned  above.  We  begin  with  the  1*  X  1*  topography 
supplied  by  NCAR.  These  data  are  first  Interpolated  onto 
the  2.5*  X  2.5*  grid  on  which  the  other  surface  fields  are 
available.  The  2.5*  X  2.5*  values  are  then  Interpolated 
onto  the  model  grid  according  to  the  method  of  Brenner 
(1987).  The  most  not'ceable  differences  are  the  higher 
elevations  over  South  America  and  around  the  Mediterranean 
basin  and  the  steeper  slopes  In  the  Himalayas. 

The  land  sea  mask  Is  a  flag  which  tells  the  model  If 
the  surface  Is  land  or  sea  for  the  purposes  of 
surface  flux  calculations  In  the  planetary  boundary  layer 
scheme.  In  the  AFGL  baseline  model,  the  land  sea  mask  Is 
based  on  the  sea  surface  temperature  field  so  that  any 
given  grid  point  can  be  either  all  land  or  all  sea.  For  this 
study  a  land  sea  mask  was  derived  which  allows  for  partial 
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land  coverage  In  a  grid  box.  The  2.5'*  X  2.5*  values  were 
derived  In  conjuctlon  with  the  topography  field  based  on  the 
original  1*  X  1“  topgraphy  data.  The  required  values  on  the 
model  grid  are  obtained  via  the  same 
Interpolation/analysis  method  used  for  topography.  The  land 
sea  mask  and  the  topography  on  the  model  grid  are  consistent 
in  the  sense  that  any  grid  point  defined  as  all  sea  cannot 
have  an  elevation  different  from  zero  before  the  spectral 
decomposition  of  the  topography  field. 

The  sea  surface  temperatures  used  in  this  study  are 
the  climatological  monthly  mean  values  which  are  available  on 
the  2.5*  X  2.5*  grid.  They  are  Interpolated  onto  the  model 
grid  using  the  same  method  as  for  topography  and  the  land 
sea  mask.  The  temperature  of  sea  Ice  is  set  to  271.15  K. 
Surface  temperature  over  land  la  a  time  dependent  field  which 
Is  obtained  by  extrapolating  downward  according  to  the  lapse 
rate  in  the  two  lowest  model  layers  (Bourke,  et  al.,  1977). 

Finally,  the  basic  values  for  the  neutral  drag 
coefficients  are  geographically  varying  (depending  upon  the 
roughness  length)  and  are  available  on  the  2.5*  X  2.5* 
grid.  These  are  also  Interpolated  onto  the  model  grid  in 
the  manner  outlined  above. 


2.2  Planetary  boundary  layer  scheme 


As  mentioned  above,  the  planetary  boundary  layer  (PBL) 
scheme  In  the  AFGL  baseline  model  consists  of  bulk  formulae 
for  surface  drag,  sensible  heat  flux,  and  evaporation,  the 
latter  two  being  allowed  only  over  the  oceans.  The  drag 
coefficient  Is  set  to  its  neutral  value  while  the 
coefficients  for  sensible  and  latent  heat  fluxes  Include  an 
additional  factor  to  account  for  wind  speed  dependence. 
Other  than  the  surface  fluxes  Into  the  lowest  layer 
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of  the  model,  there  la  no  vertical  diffusion  allowed.  Since 
our  major  interest  In  this  research  Is  air  sea  Interaction, 
we  felt  It  necessary  to  Improve  upon  the  baseline  PEL 
scheme.  However,  due  to  computer  resource  limitations,  we 
were  not  able  to  adapt  an  overly  complex  scheme  such  as  the 
one  developed  by  Mahrt  et  al.  11984)  for  AFGL,  Instead  we 
settled  on  the  scheme  described  below  which  we  feel  Is  a 
reasonable  compromise  between  being  state  of  the  art  yet 
computationally  efficient. 


Vertical  diffusion  Is  Incorporated  In  the  model  following 
the  linearized  formulation  of  Bourke  et  al.  (1977).  A  term  of 
the  form 

e 

p.  ^ 


where  g 

la  gravity. 

p.  Is  the 

surface  pressure. 

and 

F  =  <T,H, 

E)  Is  added 

to  the  momentum. 

thermodynamic, 

and 

moisture 

equations. 

respectively. 

For 

a  1  1  1  eve  1 s  except 

the 

surface. 

these  fluxes 

are  given  by 

T  = 

■5V 

Q  K,  - 

5  e 

( 1 ) 

H  = 

■je 

p  Kv  - 

3<r 

(2) 

E  = 

'^Q 

p  K,  — - 

(3) 

where  V  is  the  horizontal  wind  vector,  6  is  the  potential 
temperature,  Q  Is  the  specific  humidity,  and  the  diffusion 
coefficient,  K, ,  la  given  by 


K,  =  pz 
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The  mixing  length,  ji,  decreases  linearly  from  the  surface  up 
to  (T  =  0.5  according  to  60((r-.5)  m  and  Is  set  to  zero  above 
or  =  0.5. 

The  surface  fluxes  are  given  by  bulk  formulae  with 
stability  dependent  coefficients.  The  respective  expressions 
are 


T  = 

-  e 

c« 

V  V. 

(5) 

H  = 

p 

Ch 

V  (ft. -Qk) 

(6) 

E  = 

p 

Co 

V  (6. -e«) 

(7) 

where  C« , 

Ch 

Co  are  the  bulk  coefficients, 

V 

magnitude  of  the  wind  speed  In  the  lowest  model  layer,  and 
the  subscripts  »  and  K  Indicate  values  at  the  surface  and  at 
the  lowest  model  layer  respectively.  The  coefficients  are 
functions  of  the  bulk  Richardson  number,  Rl,  of  the  lowest 
half  layer  which  is  defined  by 
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In  determining  the 

value  of 

h  it  was 

found 

that 

the 

Arakawa  vertical 

differencing  scheme  produces  rather  large  errors  due  to  the 
use  of  an  integral  constraint  (Brenner,  1988).  For  this 
reason,  h  in  the  PBL  scheme  is  computed  according  to  the 
discrete  hydrostatic  equation  described  by  Gordon  and  Stern 
(1982).  For  unstable  stratification  (Ri  <  0),  the  bulk 
coefficients  are  given  by 


7 
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Eqs.  (9)  Is  taken  from  Deardorff  (1968)  and  the 
constant,  a,  is  given  by  a  =  0.83  Co « ' “ ■ ‘ . 
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the  neutral  drag  coefficient,  Cd n ,  Is  taken  from  the  fixed 
field  values  supplied  on  the  2.5*  X  2.5*  grid  Interpolated  on 
to  the  model  grid.  Over  water,  the  values  are  given  as  a 
function  of  the  surface  wind  speed 

Com  =  CO. 43  +  0.097  V]  X  10-^  (11) 

as  In  Geernaert  et  al.  (1966). 

For  stable  stratification  (R1  1  0),  the  bulk  coefficients 
are  given  by 


C,  =  Com  exp  (-R1)  (12) 

1.35  Cm 

Cm  =  Ct  = -  (13). 


1.  +1.93  R1 

Eq.  (12)  Is  taken  from  Mahrt  (1987)  and  eq.  (13)  Is  from 
Hansen  et  al.  (1983). 

Finally,  the  surface  wind  Is  taken  as  the  wind  of  the 
lowest  model  layer  rotated  towards  lower  pressure.  The 
turning  angle,  a,  depends  upon  the  stability  and  type  of 
surface.  Over  land  the  values  are  30*  for  Ri  2  0  and  10*  for 
R1  <  0.  Over  the  sea  the  values  are  20*  and  0*  for  stable  and 
unstable  stratification  respectively  ( 0 t to-B 1 e 1 sner ,  et  al. 
1982)  . 

2.3  Holst  convective  scheme 


The  effects  of  moist  convection  are  Included  In  the 
AFGL  baseline  model  according  to  the  cumulus  paramter 1 zat 1  on 
scheme  of  Kuo  (1965).  In  this  scheme,  a  cloud  Is  formed 
whenever  the  atmosphere  Is  unstable  and  the  vertically 
Integrated  horizontal  moisture  covergence  exceeds  a  certain 
threshold  value. 


Recently,  Soong  et  al.  (1865)  added  several  Improvements 
to  the  Kuo  scheme  in  support  of  the  global  modelling  research 
at  AFGL.  These  Improvements  Include  the  predlcltlon  of  the 
height  of  the  cloud  tops  and  the  incorporation  of  the  effects 
of  entrainment  on  the  cloud  temperature  and  specific  humidity 
profiles.  For  our  study,  we  have  adapted  this  modified 
cumulus  parameterization  scheme. 


2.4  Diagnostic  cloud  scheme 


In  the  present  version  of  the  model  there  is  no 
radiative  transfer  routine  so  that  the  main  purpose  of  the 
cloud  scheme  is  to  serve  as  a  convenient  way  of  representing 
the  humidity  field  at  any  particular  forecast  time.  Thus  we 
are  concerned  only  with  specifying  cloud  cover  but  not  the 
radiative  properties  of  the  clouds.  Ue  are  also  not 
concerned  here  with  the  complex  problem  of  c 1 oud- rad lat 1  on 
interaction.  For  this  work  we  use  the  cloud  scheme  devised 
by  Ge I eyn  et  al.  (1982)  In  which  the  cloud  amount,  C,  in  any 
model  layer  Is  determined  as  a  function  of  the  relative 
humidity  and  is  given  by 

RH-RH. 

C  =  (max  ( - ,  0.0))2  (14) 

1. -RHc 

where  RH  is  the  layer  relative  humidity  and  RH^  is  a 
critical  value  which  is  a  function  of  the  vertical  coordinate 
given  by 

RH.  =  1  -  2a(l-(r)tl  +  J3(ij-ii)]  (15). 

The  total  cloud  cover  In  a  column  Is  specified  by  maximum 
overlap,  l.e.  the  total  is  set  equal  to  the  maximum  layer 
value  in  the  column. 
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3.  EXPERIMENTS  AND  RESULTS 


In  this  section  we  describe  the  experiments  that  ue 
conducted  and  their  results.  In  view  of  our  limited  computer 
resources  we  found  It  necessary  to  choose  between  the 
number  of  experiments  and  the  spatial  resolution  of  the 
model.  In  order  to  obtain  meaningful  statistics  for  the 
significance  of  our  results,  it  was  preferable  to  run  as 
many  test  forecasts  as  possible.  Thus  for  most  of  our 
work  we  used  a  resolution  of  nine  equal  thickness  sigma 
layers  in  the  vertical  and  rhomboldal  15  truncation  in  the 
horizontal.  The  corresponding  transform  grid  is  40  latitude 
by  48  longitude  points  <4.4”  X  7.5').  This  is  comparable  to 
the  resolution  used  in  many  spectral  climate  models  (e.g. 
Pitcher  et  al.,  1983).  We  will  refer  to  these  cases  as  low 
resolution  and  designate  them  as  15R9.  All  of  the  experiments 
presented  below  were  run  on  either  a  Digital  VAX  11/750  or  a 
Mlcrovax  II.  With  a  time  step  of  45  min,  the  15R9  model 
requires  62  min  of  CPU  time  per  forecast  day.  In  one 
particular  case,  we  Increased  the  horizontal  resolution  to 
rhomboldal  23  corresponding  to  a  transform  grid  of  60 
latitudes  by  72  longitudes  (3°  X  5’).  This  will  be  designated 
23R9. 

3. 1  Data  sources 


For  the  experiments  presented  below,  we  obtained  data 
from  the  European  Center  for  Medium  Range  Weather  Forecasts 
(ECMWF),  An  extensive  archive  of  dally  global  analyses 


(OOZ  and  12Z)  beginning  from 
ECMWF.  These  data  sets  consist 
wind  components,  geopotential 


1980  Is  maintained  at 
of  initialized  values  of 
height,  temperature. 
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relative  humidity,  and  vertical  velocity.  The  data  are 
available  on  a  2.5*  X  2.5*  grid  at  the  following  seven 
mandatory  pressure  levels;  1000,  850,  700,  500,  300,  200, 
and  100  mb. 

The  data  can  be  requested  in  a  choice  of  tape 
formats.  Ue  orignally  ordered  a  tape  written  in  a 
compressed  binary  mode  in  which  up  to  three  months  of 
data  (either  OOZ  or  122)  can  be  stored  on  a  single  nine 
track  tape  at  1600  bpl.  After  several  frustrating  and 
unsuccessful  attempts  to  read  this  tape  we  reodered  the 
data  in  the  UMQ-GRID  format  which  is  a  compressed  code 
written  in  ASCII  characters.  In  this  format,  up  to  23  days 
of  data  can  be  stored  on  a  single  tape.  Reading  the  tape 
was  fairly  straightforward,  but  it  took  us  several  weeks  to 
write  the  appropriate  software  to  remove  extraneous  control 
characters  and  to  decode  the  data.  In  the  end,  ue  obtained 
the  122  data  for  the  forty  five  day  period  from  17  October 
through  30  November  1986.  Ue  chose  this  period  since  it 
Included  several  storms  which  crossed  the  Mediterranean 
and  brought  significant  rain  to  Israel  and  the  surrounding 
area . 


3.2  Low  resolution  experiments 


In  this  section  we  present  the  results  of  a  series  of  low 
resolution  (15R9)  experiments  in  which  we  examined  the 
sensitivity  of  the  model  produced  precipitation  to  changes  In 
the  sea  surface  temperatures  of  the  eastern  Mediterranean.  We 
ran  a  series  of  eleven  96  hour  forecasts  covering  successive 
four  day  periods  beginning  from  122  on  17  October  1986. 
Thus,  the  Initial  data  for  the  eleven  forecasts  consisted  of 
the  122  analyses  on  17,  21,  25,  29  October,  2,  6,  10,  14,  18, 
22,  and  26  November  1986.  From  each  of  these  Initial 
times, we  ran  two  forecasts.  In  the  control  forecast 
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(designated  CONT)  we  used  climatological  monthly  mean  sea 
surface  temperatures.  In  the  anomaly  forecast  (designated 
ANOM)  we  Increased  the  sea  surface  temperature  of  the  eastern 
Mediterranean  by  2*C  In  the  Ionian  Sea  and  by  2.5*C  In  the 
Levantine  Sea.  The  value  of  2.5*C  Is  four  times  the  standard 
deviation  of  the  multiyear  observed  winter  surface 
temperature  of  the  eastern  Levantine  Sea  (Hecht  et  al., 
1988).  Thus  In  our  experiments,  this  anomaly  Is  meant  to 
represent  an  upper  limit  on  the  range  of  the  interannual 
variability  of  the  eastern  Mediterranean  (l.e.  ±2  standard 
deviations).  Ue  chose  a  forecast  period  of  96  hours  for  two 
reasons.  First,  this  is  close  to  the  skillful  forecast  limit 
of  this  particular  model  configuration  (l.e,  low  resolution 
with  no  radiative  transfer  routine).  Second,  we  felt  that 
this  was  a  long  enough  period  to  allow  for  the  spin  up  of  the 
divergence  and  moisture  fields  In  the  model. 

Our  main  Interest  In  this  study  Is  the  sensitivity  of  the 
model  produced  precipitation  to  changes  in  sea  surface 
temperatures.  Thus  we  will  only  briefly  address  the  Issue  of 
performance  and  verification  In  order  to  demonstrate  the 
capabilities  of  the  model.  To  this  end,  In  Figure  1  we  show 
the  global  root  mean  square  (RMS)  1000-200  mb  height  error 
for  the  eleven  control  forecasts.  If  we  use  persistence  as  a 
measure  of  the  performance  limit,  then  at  four  days  the 
model  still  shows  some  skill.  Extrapolation  of  these  two 
curves  Indicates  that  the  model  will  no  longer  be  skillful 
(l.e.  It  would  lose  to  persistence)  at  about  4.5  days.  With 
this  as  background,  we  now  proceed  to  the  main  results  of 
this  study,  l.e.  the  precipitation  and  cloud  forecasts.  We 
again  emphasize  that  our  Interest  is  the  sensitivity  of  the 
model  and  so  we  will  focus  our  discussion  on  the  model 
results  and  not  on  verification  against  observed  values. 

Figure  2  shows  the  twelve  day  (three  forecasts)  running 
mean  rainfall  In  mm  day*  for  the  region  bounded  by 
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Figure  1:  Root  mean  square  (RMS)  height  error  for 
the  layer  1000-200  mb  for  the  eleven  control 
forecasts  with  the  1SR9  model. 


Figure  2:  Twelve  day  running  mean  rainfall  Ir 
mm  day-*  over  the  region  20*-50*K,  20*-45*E. 
Solid  line  -  control  forecasts;  dashed  line  - 
anomaly  forecasts. 


20*-50*N,  20*-45°E.  For  each  of  the  forty  four  forecast 
days,  the  anomaly  rainfall  exceeds  the  corresponding  control 
amount.  The  pronounced  peaks  on  22  October,  2  November,  and  8 
November  correspond  to  three  major  storms  that  passed  through 
the  region.  From  9-24  November  there  were  no  storms  as 
Indicated  by  the  sharp  decrease  In  the  mean  rainfall 
during  that  period.  The  trend  reverses  on  26  November  when 
another  storm  appeared  in  the  region.  The  largest  differences 
between  the  control  and  anomaly  amounts  tend  to  be  at  the 
times  of  the  most  Intense  rainfall  (l.e.  at  the  peaks  In  the 
curves ) . 

In  Figure  3  we  show  the  total  rainfall  in  mm  day*  for 
the  entire  forty  four  days  of  the  eleven  control  forecasts. 
There  are  three  regions  where  the  precipitation  field 
exhibits  relative  maxima.  The  most  intense  is  located  In  the 
northwest  corner  of  the  map  off  the  Atlantic  coast  of  France. 
This  is  associated  with  North  Atlantic  storms  which  cross 
central  and  northern  Europe.  In  general,  these  storms  do  not 
affect  the  eastern  Mediterranean  region.  The  second 
maximum  appears  as  a  tongue  of  precipitation  greater  than 
3  mm  day*  over  the  Adriatic  Sea  and  Yugoslavia.  The  third 
maximum  occurs  over  the  eastern  end  of  the  Mediterranean 
(values  greater  than  2  mm  day*  )  and  extends  eastward  Into 
Iraq  (off  the  edge  of  the  map).  The  regions  of  minimum 
rainfall  Include  North  Africa  and  the  part  of  the  Soviet 
Union  Just  north  of  the  Black  Sea  as  Indicated  by  the  closed 
1  mm  day*  contour. 

The  changes  In  the  precipitation  field  for  the  entire 
forty  four  day  period  Is  shown  In  Figure  4  where  we  have 
plotted  the  values  of  ANOM-CONT.  The  shading  Indicates 
statistical  significance  at  or  above  the  90)(  level  based  on  a 
t-test  for  paired  comparisons.  The  significant  increase  In 
precipitation  Is  confined  mainly  to  the  Levantine  Sea  and  the 
surrounding  area.  Maximum  values  of  0.6  mm  day*  or  more 
appear  as  a  band  which  follows  the  coast  of  Turkey.  Over 
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Figure  3:  Total  rainfall  In  ram  day'*  for  the 
entire  forty  four  day  period  of  the  eleven  control 
forecasts.  The  contour  Interval  Is  1  ram  day'  . 


Figure  4:  ANOM-CONT  rainfall  difference  for  the 
same  period  as  in  Figure  3.  Shading  Indicates 
statistical  significance  at  or  above  the  90X  level. 
The  contour  interval  Is  0.2  mm  day'. 
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central  and  northern  Israel  and  Lebanon,  the  Increase  of  0.4 
mm  day*  corresponds  to  15-20*  of  the  total  rainfall  for 
that  area.  There  Is  a  slight  decrease  In  the  rainfall  across 
eastern  Europe,  over  Algeria,  and  over  parts  of  Spain  and 
France.  However  these  decreases  are  not  significant  except 
for  the  two  very  small  areas  Indicated  on  the  map. 

The  obvious  source  of  moisture  for  storms  that  enter  into 
or  develop  over  the  Mediterranean  region  must  be  the  sea. 
The  role  of  the  Mediterranean  as  a  source  of  evaporation  is 
further  enhanced  by  its  nature  of  being  a  warm  semi -enc 1 osed 
sea.  This  effect  is  analogous  to  the  role  of  the  Great  Lakes 
in  enhancing  winter  snowfall  over  the  so  called  "snow  belt" 
as  cold  storms  from  Canada  cross  over  the  relatively  warm 
lakes.  In  Figure  5  we  show  the  total  latent  heat  flux  from 
the  sea  (expressed  as  evaporation  in  mm  day*)  for  the 
eleven  control  forecasts.  The  eastern  Mediterranean  with  its 
relatively  high  surface  temperatures  stands  out  with  a 

maximum  of  over  6  mm  day*  to  the  south  of  Cyprus.  The 
western  Mediterranean  appears  to  be  of  much  less  importance. 
The  values  for  evaporation  are  reasonable  compared  to  those 
computed  by  Bunker  (1972)  based  on  ship  observations.  The 

values  of  the  model  sensible  heat  flux  from  the  sea  (not 

shown  here)  are  typically  less  than  20*  of  the  latent  heat 
flux.  Once  again  this  seems  reasonable  compared  to  the 

computations  of  Bunker  (1972). 

In  Figure  6  we  show  the  increase  in  evaporation  for  the 
eleven  anomaly  forecasts  in  mm  day*.  It  is  not  surprising 
to  find  that  the  maximum  Increase  occurs  over  the  eastern 
Mediterranean.  The  percent  increase  in  evaporation  is 
comparable  to  the  percent  Increase  in  rainfall  over  the 
region  thus  indicating  that  most  of  the  additional  moisture 
in  the  anomaly  cases  fails  out  as  precipitation  locally. 
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Figure  6:  Increase  In  evaporation  for  the  eleven 
anomaly  forecasts  (ANOM-CONT).  The  contour  Interval 
Is  0.4  mm  day  ‘  . 


In  order  to  further  understand  the  physical  processes 
responsible  for  this  Increase  In  precipitation  ue  compared 
the  mean  circulation  for  the  control  and  anomaly  runs  In  the 
lowest  model  layer.  We  found  that  the  wind  components  and  the 
height  field  did  not  differ  substantially  (e.g.  ANOM-CONT 
differences  In  the  wind  components  were  at  most  a  few  tenths 
of  a  meter  per  second).  However,  the  temperature  of  the 
lowest  model  layer  did  show  significant  changes  as  Indicated 
In  Figure  7.  Over  the  area  of  the  Imposed  sea  surface 
temperature  anomaly,  the  mean  temperature  for  the  eleven 
anomaly  forecasts  was  at  least  0.5  K  warmer  than  In  the 
corresponding  control  cases.  This  Is  due  to  the  change  In  the 
sensible  heat  flux  from  the  warmer  sea  in  the  anomaly  runs. 
The  maximum  sensible  heat  flux  to  the  south  of  Cyprus 
Increased  from  44  W  m"  ^  for  the  control  runs  to  67  U  m" ^ 
in  the  anomaly  runs. 


Figure  7:  ANOM-CONT  difference  In  the  temperature 
of  the  lowest  model  layer  for  the  entire  forty 
four  day  period.  The  contour  interval  Is  0.2  K. 
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Figure  8:  Large  scale  rainfall  for  the  entire  forty 
four  day  period  of  the  eleven  control  runs.  The 
contour  Interval  Is  1  mm  day*. 


Figure  9:  As  In  Figure  6  except  for  the  convective 
rainfall.  The  contour  Interval  Is  0.5  mm  day*. 
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This  Increase  In  the  temperature  of  the  lowest  layer  led 
us  to  suspect  that  perhaps  the  static  stability  decreased 
and  as  a  result,  the  convective  activity  increased  in  the 
anomaly  runs.  Ue  were  able  to  check  this  since  the  model 
computes  and  saves  the  large  scale  and  convective  rainfall 
separately.  Figures  8  and  9  show  the  large  scale  and  the 
convective  components,  respectively,  of  the  total  control 
rainfall  shown  In  Figure  3.  Recalling  the  three  precipitation 
maxima  from  Figure  3,  we  can  now  see  that  the  maximum  in  the 
northwest  corner  of  the  map  is  due  mainly  to  large  scale 
precipitation.  The  land  based  part  of  the  high  rainfall 
tongue  that  appeared  over  Yugoslavia  Is  mainly  due  to  the 
large  scale  while  the  part  of  this  tongue  over  the  Adriatic 
Sea  Is  due  to  convective  actltlvty.  It  appears  that  the 
former  Is  related  to  orographic  lifting  while  the  latter  is 
related  to  convective  activity  over  the  relatively  warm  sea. 
The  third  maximum  which  extended  from  the  Levantine  Sea 
across  Israel,  Lebanon,  Syria  and  Into  Iraq  is  due  almost 
exclusively  to  convective  rainfall.  There  is  also  one 
additional  maximum  of  covective  activity  over  Algeria  and 
Libya  which  did  not  stand  out  In  the  total  rainfall  map  In 
Figure  3.  One  final  point  to  note  In  Figure  9  Is  that  the 
convective  rainfall  occurs  In  a  band  that  Is  oriented  along 
the  major  axis  of  the  eastern  Mediterranean.  It  seems  that 
the  convective  activity  Is  triggered  by  the  temperature 
contrast  between  the  relatively  warm  sea  and  the  cold  land 
surfaces  to  the  north. 

We  conclude  this  section  with  Figures  10  and  11  which  show 
the  Increase  In  large  scale  and  convective  rainfall, 
respectively,  for  the  anomaly  forecasts.  These  are  the  two 
components  of  the  total  change  that  was  shown  In  Figure  4.  It 
Is  not  surprising  to  note  that  the  major  changes  appear  over 
the  eastern  Mediterranean  In  both  components  of  the  rainfall. 
This  Is  due  to  the  increase  In  available  moisture  by  the 
enhanced  evaporation  In  the  anomaly  forecasts.  It  is 
interesting  to  note  that  the  convective  rainfall  contributes 
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Figure  10:  ANOM-CONT  difference  In  large  scale 
rainfall  for  the  eleven  cases.  The  contour  interval 
is  0.2  mm  day- ‘ . 


Figure  ll!  As  In  Figure  10  except  for  convective 
ra Inf  a  1 1 . 


two  thirds  of  the  total  Increase  while  the  large  scale 
accounts  for  only  one  third  of  the  total  Increase.  Thus  we 
can  see  that  the  Increase  In  rainfall  over  this  region  caused 
by  a  warm  sea  surface  temperature  anomaly  in  the  eastern 
Mediterranean  Is  due  to  Increased  moisture  availability 
through  enhanced  evaporation  and  Increased  convective 
activity  triggered  by  heating  of  the  atmosphere  from  below  by 
an  enhanced  sensible  heat  flux. 


3.3  Case  study  -  6  Novenber  1986 


In  order  to  further  examine  the  sensitivity  of  the  model 
we  chose  one  case  from  the  eleven  discussed  above  and  ran  a 
series  of  forecasts  with  a  range  of  values  for  the  sea 
surface  temperature  anomaly.  For  this  case,  we  also  ran 
several  forecasts  at  a  higher  horizontal  resolution  (23R9) 


Figure  12;  Global  RMS  1000-200  mb  height  error  for 
the  control  forecasts  from  12Z  6  November  1986. 
Solid  line  -  15R9,  dotted  line  -  23R9i  and  dashed 
line  -  persistence. 
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for  comparison  with  the  low  resolution  runs.  The  forecasts 
for  this  case  study  were  run  from  the  Initial  data  at  12Z  on 
6  November  1S86.  This  particular  case  was  chosen  since  it 
corresponds  to  the  highest  rainfall  rate  In  Figure  2.  As  In 
the  previous  section,  before  presenting  the  main  results,  we 
briefly  address  the  performance  of  the  model.  In  Figure  12  we 
show  the  global  RMS  1000-200  height  error  for  the  15R9  and 
23R9  control  forecasts  for  the  6  November  case.  Ue  also 
include  the  persistence  curve  for  comparison.  As  expected, 
the  23R9  error  is  smaller  than  the  15R9  error.  Extrapolation 
of  these  curves  indicates  that  the  higher  resolution  forecast 
retains  Its  skill  for  about  one  day  longer  than  the  lower 
resolution  forecast. 

The  synoptic  situation  on  8  November  (day  of  maximum 
predicted  rainfall)  Is  depicted  In  Figures  13  and  14  where  ue 
show  the  48  hour  23R9  forecast  1000  mb  and  500  mb  height 
fields,  respectively.  This  Is  a  very  typical  winter  rain 
producing  situation  for  the  eastern  Mediterranean  and  the 
surrounding  area.  The  main  surface  low  Is  centered  over 
eastern  Turkey  (off  the  edge  of  the  map)  with  a  cold  front 
extending  to  the  southwest  across  the  Levantine  Sea.  There  Is 
also  a  hint  of  a  secondary  Cyprus  low  developing.  Strong 
upper  level  support  for  this  system  Is  apparent  In  the  500  mb 
map.  The  trough  extends  across  the  Levantine  Sea  and  the 
strong  northerly  winds  supply  this  system  with  cold  polar  air 
which  Is  a  major  ingredient  necessary  for  heavy  rain  In  this 
region. 

In  Figures  15  and  16  we  show  the  total  rainfall  for  the 
15R9  and  23R9  control  forecasts  from  6  November.  The  general 
pattern  Is  similar  for  both  forecasts  although  there  are 
certainly  some  noticeable  differences  In  the  details.  Both 
resolutions  correctly  show  the  heaviest  rainfall  associated 
with  the  cold  front  over  Turkey  and  the  eastern  Mediterranean 
with  maximum  values  of  over  6  mm  day* ■ .  However,  In  the 
23R9  run  the  sharpness  of  the  cold  front  is  better  defined 
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Figure  13:  1000  mb  height  48  hour  CONT  forecast  for 
8  November  1986  with  the  23R9  model.  The  contour 
Interval  Is  2  decameters. 
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Figure  14:  As  In  figure  13  except  500  mb  heights. 
The  contour  Interval  Is  6  decameters. 
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Figure  16:  As  In  Figure  15  except  for  the  23R9  CONT 
forecast. 


and  the  heavy  rainfall  (>  4  mm  day*)  Is  restricted  to  a 
much  narrower  band  around  the  front.  The  other  noticeable 
differences  are  the  appearance  of  two  regions  of  higher 
rainfall  along  the  no. them  edge  of  the  map  on  the  23R9 
forecast,  and  the  region  of  no  rainfall  centered  over  Italy 
Is  much  broader  in  the  23R9  map. 

Figures  17  and  18  show  the  changes  In  rainfall  for  the 
15R9  and  23R9  anomaly  forecasts,  respectively.  As  In  the  case 
of  total  rainfall,  the  overall  patterns  for  the  two 
resolutions  are  similar  although  they  do  differ  in  detail. 
Both  exhibit  an  Increase  confined  mainly  to  the  area  of  the 
Imposed  sea  surface  temperature  anomaly  and  extending 
eastward.  Thus  In  the  anomaly  runs,  the  area  of  maximum 
rainfall  Is  shifted  southward  and  eastward.  The  maximum 

increase  Is  larger  by  0.4  mm  day*  In  the  15R9  case.  This 
difference,  however.  Is  related  to  the  effective  strengh  of 
the  sea  surface  temperature  anomaly  which  Is  somewhat  weaker 
In  the  23R9  case.  Due  to  the  difference  between  the  size  of 
the  grid  boxes  in  the  two  resolutions  and  the  resulting 
difference  in  the  area  of  the  sea  (l.e.  the  land  sea  mask). 
It  is  difficult  to  apply  exactly  the  same  strength  anomaly  to 
the  two  models.  The  other  difference  between  the  two  figures 
is  the  region  of  decreased  rainfall  (indicated  by  shading). 
In  the  15R9  case,  it  appears  as  a  broad  T-shaped  feature 
while  in  the  23R9  case  it  broken  up  into  three  separate 
areas.  Nevertheless,  the  patterns  do  resemble  one  another. 
Although  It  Is  difficult  to  generalize  on  the  basis  of  one 
test  case,  we  feel  that  the  qualitative  similarities  between 
the  low  and  higher  resolution  results  are  encouraging  and 
further  support  the  conclusions  we  reached  in  Section  3.2 
which  were  based  only  on  the  low  resolution  model. 

We  noted  above  that  the  Intensity  of  the  sea  surface 
temperature  anomaly  used  in  these  experiments  was  meant  to 
represent  an  upper  limit  on  the  observed  range  of  Interannual 
var lab  I i 1 1 ty .  The  next  question  of  course  is  how  important  is 
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V 


the  magnitude  of  the  anomaly  Is  changing  the  rainfall.  To 
address  this  question  we  ran  several  additional  experiments 
for  the  6  November  case  in  which  we  added  sea  surface 
temperature  anomalies  in  the  same  locations  but  of  different 
magnitudes.  For  the  15R9  model,  In  addition  to  the  standard 
value  of  +2.5*0  in  the  Levantine  Sea,  we  also  tried  anomalies 
of  -2.5*C,  +1.0*C.  and  +4.0*.  For  the  23R9  model  we  also 
tried  a  value  of  -2.5*C.  The  results  are  shown  in  Figure  19 


Figure  19:  Total  rainfall  over  the  area  20*-50*N, 
20”-45*E  for  the  6  November  case  as  a  function  of 
the  magnitude  of  the  sea  surface  anomaly  in  the 
Levantine  Sea. 

where  we  show  the  total  rainfall  over  the  region  20”-50*N, 
20*-45*E  as  a  function  of  the  magnitude  of  the  sea  surface 
temperature  anomaly  in  the  Levantine  Sea.  From  these  two 
curves  we  can  clearly  see  the  general  trend  for  the  rainfall 
to  Increase  with  the  magnitude  of  the  sea  surface  temperature 
anomaly.  For  the  15R9  case,  the  fluctuations  of  the  area 
averaged  rainfall  are  about  ±15X  of  the  CONT  (i.e.  zero 
anomaly)  case.  At  the  upper  end  of  the  curve,  there  is  a  hint 
that  the  rate  of  increase  will  gradually  level  off  and 
possibly  reach  some  asymptotic  value.  For  the  23R9  runs,  the 
fluctuations  are  much  less  dramatic.  Furthermore,  the 
consistently  lower  values  in  the  23R9  runs,  as  compared  to 


-27- 


1BR9,  are  related  to  the  differences  In  resolution  and  the 
ability  of  the  23R9  model  to  better  represent  smaller  scale 
variations  in  the  humidity  field.  In  any  event,  these  results 
indicate  that  influence  of  the  eastern  Mediterranean  on  the 
rainfall  is  very  localized  and  is  confined  mainly  to  the  area 
over  and  immediately  surrounding  the  sea. 

In  Figure  20  we  show  the  change  in  rainfall  for  the  -2.5*0 
anomaly  run  with  the  23R9  model.  As  in  the  other  cases,  the 
maximum  change  (decrease  in  this  case)  occurs  over  the  area 
of  the  Imposed  sea  surface  temperature  anomaly.  The  magnitude 
of  the  decrease  here  is  not  as  large  as  the  corresponding 
increase  for  the  +2.5*0  anomaly  (compare  to  Figure  18).  Ue 
note  that  the  cold  anomaly  is  quite  effective  in  surpressing 
most  of  the  rainfall  over  southern  Israel,  Jordan  and  Saudla 
Arabia.  For  this  particular  case,  the  cold  anomaly  has  caused 
the  band  of  maximum  precipitation  to  become  narrower  and 
shifted  to  the  north. 


Figure  20:  ANOM-OONT  rainfall  difference  for  the 
23R9  cold  (-2.5*0)  anomaly  forecast  from  6  November 
1966.  Shaded  areas  are  negative.  The  contour 
interval  is  1  mm  day-'. 
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We  conclude  this  section  with  some  results  of  the  cloud 
forecasts  for  the  23R9  model.  Figure  21  shows  the  48  hour 
control  forecast  for  total  cloud  cover  at  12Z  on  8  November 
1986.  A  comparison  with  the  synoptic  situation  as  shown  In 
Figures  13  and  14  indicates  a  reasonable  distribution  of 
cloud  cover.  The  two  bands  of  heavy  clouds  -  along  the 
eastern  edge  and  In  the  northwest  quadrant  of  the  map  -  are 
associated  with  the  two  lows  that  appear  In  the  1000  mb 
height  forecast  (Figure  13).  The  region  of  clear  skies 
associated  with  the  high  over  northwest  Africa  is  also  well 
predicted.  The  cloud  cover  for  the  anomaly  run  is  very 
similar  and  is  not  shown  here.  In  Figure  22,  we  show  the 
cloud  forecast  for  the  cold  anomaly  (-2.5*0  run  which 
differs  substantially  from  the  control  forecast  over  the 
eastern  Mediterranean  region.  The  maximum  cloud  cover  has 
been  reduced  and  the  band  of  heavy  clouds  has  shifted  to  the 
northwest.  The  southern  half  of  the  Levantine  Sea  Is  now 
clear  as  compared  to  partly  cloud  or  cloudy  In  the  control 
f  oreoas  t . 

4.  CONCLUSIONS 


We  have  presented  the  results  of  a  series  of  numerical 
experiments  designed  to  study  the  Influence  of  the  eastern 
Mediterranean  Sea  on  the  precipitation  In  the  Middle  East.  We 
concentrated  our  attention  on  the  period  from  17  October 
through  30  November  1986  which  was  a  relatively  wet  period 
for  this  region. 

In  the  first  set  of  experiments  we  ran  a  series  of 
consecutive  96  hour  forecasts  with  a  low  resolution  (15R9) 
modified  version  of  the  AFGL  global  model.  In  one  set  of  runs 
we  used  climatological  monthly  mean  sea  suface  temperatures 
while  In  the  second  set  we  introduced  a  warm  sea  surface 
temperature  anomaly  of  2*-2.5*C  In  the  eastern  Mediterranean. 
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In  the  control  forecasts,  It  was  found  that  the  rainfall  over 
the  eastern  Mediterranean  region  and  the  Middle  East  was 
produced  almost  exclusively  by  the  moist  convection  scheme  as 
compared  to  the  rainfall  over  western  Europe  which  was 
produced  by  the  large  scale  scheme.  The  convective  rainfall 
occurs  In  a  band  which  Is  oriented  along  the  major  axis  of 
the  eastern  Mediterranean  and  is  related  to  the  temperature 
contrast  between  the  relatively  warm  sea  and  the  cold  land 
surface  to  the  north. 

The  Introduction  of  a  warm  sea  surface  temperature  anomaly 
to  the  eastern  Mediterranean  caused  an  Increase  In  the 
rainfall  over  the  region.  The  Increase  appeared  primarily  In 
the  convective  precipitation.  This  Increase  is  due  to  the 
combined  effects  of  enhanced  sensible  and  latent  heat  fluxes 
from  the  warmer  anomaly  sea  surface.  The  latent  heat  flux 
provides  additional  moisture  while  the  sensible  heat  flux 
warms  and  destabilizes  the  lower  atmosphere  and  leads  to  an 
increase  In  convective  activity. 

Ue  then  conducted  a  case  study  for  the  forecast  from 
6  November  1986.  For  this  case  ue  ran  a  series  of  experiments 
In  which  we  varied  the  magnitude  of  the  sea  surface 
temperature  anomaly.  For  the  values  tested,  it  was  found  that 
the  rainfall  Increased  monotonlca I ly  with  the  magnitude  of 
the  anomaly.  At  the  higher  end  of  the  range  there  was  some 
indication  that  the  rainfall  may  reach  some  asymptotic  limit. 
Ue  also  ran  several  experiments  In  which  the  horizontal 
resolution  of  the  model  was  increased.  A  comparison  of  the 
high  and  low  resolution  forecasts  shows  good  qualitative 
agreement  in  the  location,  Intensity,  and  pattern  of  the 
precipitation  field.  However  there  are  some  differences  In 
the  details  which  are  due  to  the  ability  of  the  higher 
resolution  model  to  better  represent  the  smaller  scale 
variations  In  the  humidity  field.  Nevertheless,  the  general 
agreement  between  the  low  and  high  resolution  runs  Is  quite 
encouraging. 


-31- 


5 


REFERENCES 


Arakawa,  A.  and  V.R.  Lamb,  1977:  Computational  design 
of  the  basic  dynamical  processes  of  the  UCLA  general 
circulation  model.  Methods  In  Computational  Physics.  Vo  1 
17,  Academic  Press,  New  York,  174-265. 

Bunker,  A.F.,  1972:  Wintertime  Interaction  of  the 
atmosphere  with  the  Mediterranean  Sea.  J.  Phvs.  Oceanogr. .  2, 
225-238. 

Bourke,  W. ,  B.  McAveney,  K.  Purl,  and  R.  Thurllng,  1977: 
Spectral  methods  In  atmospheric  modelling.  Methods  In 
Computational  Physics.  Vol.  17,  Academic  Press,  New  York. 

Brenner,  S. ,  1987:  Derivation  of  consistent  surface  fixed 
fields  for  use  In  a  global  NWP  model.  AFOSR  Scl.  Rep.  AFOSR- 
87-0066-001,  19pp.  GL-TR-89-0U2 

- ,  1988:  A  comparison  of  various  numerical  solutions 

of  the  hydrostaclc  equation.  Reprints  of  the  Eighth 
Conference  on  Numerical  Weather  Prediction.  22-26  February 
1988.  Baltimore.  MD.  American  Meteorological  Society,  438- 
44 1‘. 


- ,  S. ,  C.H.  Yang  and  S.  Yee,  1982:  The  AFGL  spectral 

model  of  the  moist  global  atmosphere:  Documentation  of  the 
baseline  version.  AFGL  Tech.  Rep.  82-0393,  65  pp.  ADA129283 

. ,  C.H.  Yang  and  K.  Mitchell,  1984:  The  AFGL  global 

spectral  model;  Expanded  resolution  baseline  version.  AFGL 
Tech.  Rep.  64-0306,  72  pp.  ADA160370 

Cressman,  G.P.,  1959:  An  operational  objective  analysis 

system.  Mon.  Wea .  Re v . .  87,  367-374. 

Deardorff,  J.W.,  1968:  Dependence  of  air  sea  transfer 
coefficients  on  bulk  stability.  J .  Geophvs .  Rea . .  73,  2549- 
2558. 

Geleyn,  J.F.,  A.  Hense  and  H.J.  Preuss,  1962:  A 
comparison  of  model  generated  radiation  fields  with 
satellite  measurements.  Beitr.  Atmos.  Phys..  55,  253-266. 

Gordon,  C.T.  and  W.F.  Stern,  1962:  A  description  of  the 
GFDL  global  spectral  model.  Mon.  Wea.  Rev. .  110,  625-644. 

Geernaert,  G.L.,  K.B.  Katsoros,  and  K.  Richter,  1986: 
Variation  of  the  drag  coefficient  and  and  its  dependence  on 
sea  state.  J .  Geophvs.  Res. .  91,  7667-7679. 

Hansen,  J.,  G.  Russell,  D.  Rind,  P.  Stone,  A.  Lads,  S. 
Lebedeff,  R.  Ruedy,  and  L.  Travis,  1983:  Efficient  three 
dimensional  global  models  for  climate  studies.  Models  1  and 
II.  Mon.  Wea.  Rev..  Ill,  609-662. 


-32- 


Hecht,  A.,  N.  Plnardl,  and  A.R.  Robinson,  1988!  Curents, 
water  masses,  eddies,  and  Jets  In  the  eastern  Mediterranean 
Levantine  Basin.  J.  Phvs.  Qceanogr. .  10,  in  press. 

Kuo,  H.L.,  1965:  On  the  formation  and  intensification  of 
tropical  cyclones  through  latent  heat  release  by  cumulus 
convection.  Jj_  Atmos.  Scl . .  22,  40-63. 

Mahrt,  L.  ,  1987:  Grid  averaged  surface  fluxes.  Mon.  Ulea. 
Rev- ■  115,  1550-1560. 

- ,  H.  Pan,  J.  Pauraler,  and  I.  Troen,  1984:  A  boundary 

layer  parameterization  for  a  general  circulation  model. 
AFGL  Tech.  Rep.  84-0063,  179  pp.  ADA144224 

Meslnger,  F. ,  1985:  The  sigma  system  problem.  Reprints  of 
the  Seventh  Conference  on  Numerical  Weather  Prediction.  17-20 
June  1985.  Montreal.  P.Q.. Canada.  American  Meteorological 
Society,  340-347. 

0 t to-B 1 e 1 sner ,  B.U.,  G.W.  Branstator,  and  D.D.  Houghton, 
1982:  A  global  low  order  spectral  general  circulation  model. 
Part  1:  formulation  and  seasonal  climatology.  J .  Atmos .  Sci . . 
39,  929-948. 

Ozsoy,  E. ,  1981:  On  the  atmospheric  factors 
affecting  the  Levantine  Sea.  ECMWF  Tech.  Rep.  25,  29  pp. 

Pitcher,  E.J.,  R.C.  Malone,  V.  Ramanathan,  M.L. 
Blackmon,  K.  Puri,  and  W.  Bourke,  1983:  January  and  July 
simulations  with  a  spectral  general  circulation  model.  J . 
Atmos ■  Scl ■ .  40,  580-604. 

Sela,  J.G.,  1980:  Spectral  modeling  at  the 
National  Meteorological  Center.  Mon .  Wea ■  Rev . .  108,  1279- 
1292. 

Soong,  S.T.,  Y.  Ogura,  and  W.S.  Kau,  1985:  A  study  of 
cumulus  parameterization  in  a  global  circulation  model.  AFGL 
Tech.  Rep.  85-0160,  113  PP* /j)ai70137 

Tzvetkov,  E. ,  1965;  Heat  storage  In  the  eastern 
Mediterranean  and  winter  precipitation  In  Israel.  1 sr .  J . 
Earth  Scl . .  34,  102-109. 

- ,  and  G.  Assaf,  1982:  The  Mediterranean  heat 

storage  and  Israeli  precipitation.  Water  Resours.  Res..  18, 
1036-1040. 

Wallace,  J.M.,  S.  Ttbaldi,  and  A.J,  Simmons,  1983: 
Reduction  of  systematic  forecast  errors  In  the  ECMWF 
model  through  the  introduction  of  an  envelope  orography. 
Quart.  J.  Roy.  Met.  Soc. .  109,  603-717, 


-33- 


